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4^ ABSTRACT 

We present a new 400 ks Chandra X-ray observation of the merging galaxy cluster Abell 
Q 2146. This deep observation reveals detailed structure associated with the major merger event 

}^ including the Mach M = 2.3 ± 0.2 bow shock ahead of the dense, ram pressure stripped sub- 

cluster core and the first known example of an upstream shock in the ICM (M = 1 .6 ± 0. 1). 
By measuring the electron temperature profile behind each shock front, we determine the 
timescale for the electron population to thermally equilibrate with the shock-heated ions. We 
^ find that the temperature profile behind the bow shock is consistent with the timescale for 

^ Coulomb coUisional equilibration and the postshock temperature is lower than expected for 

instant shock-heating of the electrons. Although hke the Bullet cluster the electron tempera- 
tures behind the upstream shock front are hotter than expected, favouring the instant heating 
model, the uncertainty on the temperature values is greater here and there is significant sub- 
structure complicating the interpretation. We also measured the width of each shock front and 
the contact discontinuity on the leading edge of the subcluster core to investigate the suppres- 
sion of transport processes in the ICM. The upstream shock is ~ 440 kpc in length but appears 
remarkably narrow over this distance with a best-fit width of only 6^3 kpc compared with the 
■T-H mean free path of 23 ± 5 kpc. The leading edge of the subcluster core is also narrow with 

sjjl an upper limit on the width of only 2 kpc separating the cool, multiphase gas at 0.5 — 2 keV 

. ^ from the shock-heated surrounding ICM at ^ 6 keV. The strong suppression of diffusion and 

conduction across this edge suggests a magnetic draping layer may have formed around the 
^ subcluster core. The deep Chandra observation has also revealed a cool, dense plume of mate- 

rial extending ^170 kpc perpendicular to the merger axis, which is likely to be the disrupted 
remnant of the primary cluster core. This asymmetry in the cluster morphology indicates the 
merger has a non-zero impact parameter. We suggest that this also explains why the SW edge 
of the subcluster core is narrow and stable over ~ 150 kpc in length but the NE edge is broad 
and being stripped of material. 

Key words: X-rays: galaxies: clusters — galaxies: clusters: Abell 2146 — intergalactic 
medium 



1 INTRODUCTION 



Galaxy clusters are forincd through mergers of smaller subclusters 
and groups, which collide at velocities of several thousand km s~' . 
* E-mail: helen.russell@uwaterloo.ca The total kinetic energy of these mergers can reach 10^ erg, a 
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significant fraction of wiiich is dissipated by large-scale shocks 
and turbulence over the merger lifetime (for a review see [Marke^ 
|vitch & Vikhlinin|2007| l. Shocks and turbulence generated by the 
merger are also expected to amplify magnetic fields in the clus- 
ter and accelerate relativistic particles. These non-thermal phenom- 
ena have been revealed through the detection of Mpc-scale syn- 
chrotron radio halos (for recent reviews see |Feretti & Giovannini| 
[2008 : Ferrari et al. 2008) and inverse-Compton hard X-ray emis- 
sion ( Fusco-Femiano, Landi & Orlandini 2005 , Rephael T, Gruber| 
[&¥lanco,1999l but see also |Wik et al. 2009 1. The combination of 
X-ray and gravitational lensing studies of merging clusters has also 
produced compelling evidence for dark matter (eg. 'Clowe, Gon-| 
[zalez & Ma rkevitch 2004; Clowe et al. 2006. Bradac et al. 2006) 
and constraints on the dark matter self-interaction cross section (eg. 
[Randall et al"l2008,) . 

Chandra's subarcsecond angular resolution revealed sharp X- 
ray surface brightness edges in merging systems. These edges cor- 
respond to cold fronts or contact discontinuities between regions 
of gas with different entropies (Markevitch et al. 2000; Vikhlinin, 
[Markevitch & Murray|2001a[ pVIarkevitch & Vikhlinin,2007| , and 
shock fronts driven by infalling subclusters. However, whilst cold 
fronts are also found in relaxed clusters and are common in clus- 



Table 1. Details of the Chandra observations analysed in this paper. 



ter cores (eg. |Owers et al. 2009 1, there are only a few confirmed 
detections of shock fronts with a sharp density edge and an un- 
ambiguous temperature jump (the Bullet cluster, [Markevitch et al. 



2002 



et al, 



et al 



Abell 520,|Markevitch et al. 2005 , two in Abell 2146, |Russell 



20T0l Abell 754, Macario'et al.|20l 1 and Abell 2744, [Owers 



201 1[ >. These surface brightness edges are key observational 



tools for studying merging systems and provide currently the only 
method for measuring the bulk velocities of the gas in the plane of 
the sky and determining the velocity and kinematics of the merger. 

Detailed observations of shock fronts and cold fronts have 
been used to probe the relatively unknown transport processes in 
the intracluster medium (ICM). [Markevitchl ( [2006^ used a deep ob- 
servation of the bow shock in the Bullet cluster to produce the first 
measurement of the electron-ion thermal equilibration timescale in 
the ICM and determine that the timescale is likely to be shorter than 
the Coulomb coUisional timescale (eg. Fox & Loeb|1997{[Marke-[ 



[vitch|2006[[Markevitch & Vikhlinin|2007) . This exciting result sug 
gests a heating process that operates faster than Coulomb collisions 
could be operating in the magnetized ICM (eg. [Schekochihin et al.[ 
[2005[ [Schekochihin et ar][2008^ . Observations of the sharp tem- 
perature and density jumps at cold fronts suggest that thermal con- 
duction and diffusion are strongly suppressed across these edges 
(eg. Ettori & Fabian 20001. A detailed study of the cold front in 
Abell 3667 found that the width of the density jump is smaller than 
the Coulomb mean free path in the ICM ( [Vikhlinin, Markevitch &[ 
[Murray|2001a) . This edge also appears sharp and stable within a 
large sector of ±30° around the symmetry axis suggesting that hy- 
drodynamic instabilities are also suppressed. The flow of the am- 
bient ICM around the dense subcluster core will stretch initially 
tangled magnetic field lines to form a draping layer with a mag- 
netic field parallel to the front (Vikhlinin, Markevitch & Murray' 
2001b; Vikhlinin & Markevitch 2002, Asai, Fukuda & Matsumoto 
2005 , Lyutikov 2006). This magnetic draping layer could provide a 
stabilising mechanism and will suppress transport processes across 
the edge. The long, straight edges of the bullet subcluster in the 
Bullet cluster also suggest a strong suppression of turbulence by 
such a stabilising magnetic layer (Markevitch 2006 j . 



The galaxy cluster Abell 2146 {z = 0.234; [Struble & Rood] 
[T999i[BBhringer et al.|20"00) is a spectacular merging system with 
two large Mach M ~ 2 shock fronts ( [Russell et al.|2010l ). The X-ray 
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morphology suggests a recent collision, 0.1 — 0.3 Gyr ago, where 
a subcluster containing a dense cool core has passed through the 
centre of a primary cluster. The subcluster is driving a bow shock 
through the ICM and is trailed by a tail of ram pressure stripped 
material. An upstream shock is also observed to be propagating in 
the opposite direction through the outskirts of the primary cluster. 
[Canning et al.[ ( [201 1^ used the line of sight velocity difference be- 
tween the subcluster and primary cluster brightest cluster galaxies 
(BCGs) with the shock velocities to estimate that the merger axis 
is inclined at only ~ 17° to the plane of the sky. This conclusion 
is also supported by the sharp surface brightness edges of the bow 
and upstream shock fronts, which would be smeared by projection 
for a larger angle to the line of sight. Abell 2146 appears to have 
undergone a simple merger between two smaller clusters, viewed 
close to side on, and therefore has a remarkably similar structure to 
the Bullet cluster ( [Markevitch et al.|2002|[Markevitch|2006| l. 

In this paper, we present results from a deep 400 ks Chandra 
observation of Abell 2146 studying the transport processes in the 
ICM. In section|2] we discuss the Chandra data reduction, analyse 
the X-ray morphology and present maps of the ICM temperature, 
density and metallicity. In section [3] we analyse the shock fronts 
in detail and compare the postshock electron temperature profiles 
behind the bow and upstream shock fronts with models for instant 
and coUisional electron-ion equilibration. In section |4] we study 
the ram pressure stripping of the subcluster core and study the sup- 
pression of diffusion and conduction across the leading edge of the 
subcluster core. We assume Ho = 70 km s"' Mpc^i , n„, = 0.3 and 
flA = 0-7, translating to a scale of 3.7 kpc per arcsec at the redshift 
z = 0.234 of Abell 2146. All errors are IcT unless otherwise noted. 



2 CHANDRA DATA ANALYSIS 
2.1 Data reduction 

Abell 2146 was observed with the Chandra ACIS-I detector for a 
total of 377 ks split into eight separate observations between Au- 
gust and October, 2010 (Table [TJ. These new observations were 
analysed together with the archival ACIS-S observations taken in 
April, 2009 ( [Russell et al.|2010[ Table [TJ. All datasets were repro- 
cessed with CIAO 4.3 and CALDB 4.4.0 provided by the Chandra 
X-ray Center (CXC). The level 1 event files were reprocessed to 
apply the latest gain and charge transfer inefficiency correction and 
then filtered to remove photons detected with bad grades. The im- 
proved background screening provided by VFAINT mode was also 
applied. Background light curves were extracted from the level 2 
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Figure 1. Raw counts image of Abell 2146 in the 0.3-7.0 keV energy band. The image has been binned by a factor of 2. 



event files of neighbouring chips for observations on ACIS-I and 
from ACIS-Sl for observations on ACIS-S3. The background light 
curves were filtered using the LC_CLEAN scrip|^ provided by M. 
Markevitch to identify periods affected by flares. There were no 
major flares in any of the observations of Abell 2146 producing a 
final cleaned exposure of 418 ks. 

The cleaned events files were then reprojected to match the po- 
sition of the obs. ID 12245 observation. Fig.[T|shows the total com- 
bined image produced by summing images in the 0.3 — 7.0 keV en- 
ergy band extracted from each individual reprojected dataset. This 
summed image was then corrected for exposure variation by divid- 
ing by the summed 1 .5 keV monochromatic exposure maps created 
for each dataset. Point sources were identified using the CIAO al- 
gorithm WAVDETECT, visually confirmed and excluded from the 



See |http:/ /cxc.harvard.edu/contrib/maxim/acisbg/| 
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analysis using elliptical apertures with radii set to five times the 
measured PSF width ( [Freeman et al.|2002l l. 

Standard blank-sky backgrounds were extracted for each chip 
in each observation, processed identically to the events file and re- 
projected to the corresponding sky position. Each blank-sky back- 
ground was normalized to match the count rate in the 9.5 — 12 keV 
energy band in the observed dataset. This correction was less than 
10% for each dataset. The normalized blank-sky background events 
files for each chip in an observation were then split to ensure each 
had the same ratio of exposure time to the observed exposure time. 
These were then summed together to produce background events 
files which covered all chips in each observation. By matching to 
the hard X-ray background count rate, we may over or underesti- 
mate the soft component of the X-ray background. The normalized 
blank sky background datasets were tested by comparison with ob- 
served background spectra extracted from source-free regions of 
the chips. We found that the normalized blank sky backgrounds 
were a close match to the observed background over the whole en- 
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ergy band. Total blank-sky background images combining all of the 
observations were generated in a similar way to the total images as 
detailed above. 

2.2 Imaging analysis 

Fig.[2](upper left) shows an exposure-corrected image of the cluster 
X-ray emission produced by combining all of the individual Chan- 
dra observations. The cluster gas is extended along the merger axis, 
NW to SE, and the bright, dense core of the subcluster is offset from 
the centre and being stripped of its material by ram pressure in the 
collision. The X-ray morphology suggests a major merger where 
the subcluster has recently passed through the centre of the primary 
cluster. There is no obvious surface brightness peak associated with 
a primary cluster core. The primary cluster may not have originally 
had a dense core or it could have been significantly disrupted in the 
collision with the subcluster core. The subcluster is observed soon 
after core passage when it has emerged from the primary core and 
is travelling towards the SE. Fig. |2] (lower right) shows a Subaru 
R-band image of the galaxy distribution (King et al. in prep.). The 
distribution of red sequence galaxies appears to separate into two 
groups, corresponding to the locations of the subcluster and the pri- 
mary cluster, which supports the interpretation of a major merger 
( [Canning etal.|2011^ . 

The two shock fronts, reported in 'Russe ll et al.| ( (2010[ >, are 
clearly visible in the unsharp-masked image (Fig.|2] upper right) 
as surface brightness edges to the SE and NW. The SE edge corre- 
sponds to the bow shock, which has formed ahead of the subcluster 
core, and can now be traced over --^ 500 kpc in length. The NW 
surface brightness edge corresponds to the upstream shock, which 
has formed in the wake of the subcluster's passage through the pri- 
mary cluster core and is travelling in the opposite direction to the 
bow shock ( [Russell et al.|2010| . The upstream shock is visible over 
'-^ 440 kpc in length and appears to have greater curvature compared 
to the bow shock. Note that in Fig. |2] (upper right) a point source 
has been removed at the NE end of the bow shock, which has the 
effect of increasing the apparent curvature of the shock front. There 
is no obvious second surface brightness peak corresponding to the 
primary cluster core and this may have been completely disrupted 
in the collision with the subcluster. 

The deeper Chandra observations have also now revealed the 
complex structure of the cool, dense subcluster core and the ram 
pressure stripped tail. Whilst the leading edge of the core is smooth, 
narrow and roughly spherical, there is a clear difference between 
the NE and SW edges of the tail (Fig.|2] lower left). The SW edge 
appears sharp and narrow over a distance of --^ 40 arcsec. In com- 
parison, the NE edge is poorly defined, it appears broader and dis- 
rupted, suggesting that the interface with the ambient medium has 
become turbulent here and instabilities could be developing. There 
is also an extended plume of emission to the SW from the subclus- 
ter tail, ~ 45 arcsec in length, which is perpendicular to the merger 
axis. These features are discussed further in section|4] 

Fig. |3] shows the subcluster core in detail and the location of 
the BCG immediately behind the X-ray peak. The deep Chandra 
dataset confirms the detection of a hard X-ray (4 — 7 keV) point 
source, at the 3C7 level, coincident with a radio point source de- 
tected in VLA 1.4 GHz archival data (NRAO/VLA Archive Sur- 
vey) and observations with the AMI Large Array at 16 GHz ( [AMI[ 
[Consortium: Rodrfguez-Gonzalvez et al. 201 Ij l. Also, Spitzer oh- 
servations of the BCG by |Quillen et al.[ 2008| l suggest that there 
is a strong contribution to the IR emission from an AGN. This ra- 
dio and X-ray point source likely corresponds to an AGN at the 



centre of the BCG. However it was difficult to determine the flux 
as the nucleus is superimposed on a bright, dense filament of gas, 
which is detected in soft X-rays with Chandra but also in Ha and 
[N II] ^Canning et al.|20TT) . We extracted the AGN source counts 
in a region of 2 arcsec radius and subtracted the cluster emission 
using a neighbouring region to the SE from 2.5 — 5 radius, which 
lay on the dense gas filament. For the energy band 2 — 7 keV, the 
point source was detected as 210 ± 30 counts above the background 
cluster emission. Assuming a powerlaw model with photon index 
r = 2 and Galactic absorption of «h = 0.03 x 10 cm^- ( Kalberlal 
[et al.|2005) , we estimated the point source luminosity in the energy 
range 2 - 10 keV to be = 1 .5 ± 0.2 x lO*^ erg s" ' . 



2.3 Spatially resolved spectroscopy 

2.3.1 Contour binning maps 

We used spatially resolved spectroscopy techniques to produce de- 
tailed maps of the projected gas properties (Fig. |4]l. The central 
--^4x4 arcmin was divided into regions using the Contour Bin- 
ning algorithm ( Sanders 2006 ), which traces the surface brightness 
variations to generate the spatial bins. For the temperature and nor- 
malization maps, regions with a signal-to-noise ratio of 32 (~ 1000 
counts) were chosen. The length of the regions was restricted to be 
at most two and a half times their width. For each Chandra observa- 
tion, we extracted a spectrum from each of the regions, subtracted 
the background spectrum and generated appropriate responses and 
ancillary responses. The spectra were grouped to contain a mini- 
mum of 20 counts per spectral charmel and restricted to the energy 
range 0.5 — 7 keV. We summed together spectra and background 
spectra for a particular region for observations on the same chip, as 
the roll angles are also comparable in each case. The response files 
were also averaged together, weighting by the fraction of the total 
counts in each observation. Each total spectrum was fitted in XSPEC 
12 ( Arnaud 1996) with an absorbed MEKAL model. The absorption 
was fixed to the Galactic value «h = 3.0 x 10 cm" -2 ( Kalberial 
[et al.|20()5) and the redshift was fixed to 0.234. The C-statistic was 
minimised in the spectral fitting ( [Cash[[T979| ). The errors are ap- 
proximately ~ 6% in emission measure and ~ 15% in temperature. 
However, the error in temperature drops to less than 10% in the sub- 
cluster core, where the temperature falls below 2 keV and the Fe L 
line emission improves temperature diagnostics. In addition, tem- 
peratures over 10 keV are poorly constrained by the energy range 
of Chandra and the errors increase to ~ 30% in these bins. 

Fig.|4](upper left) shows the distribution of the projected emis- 
sion measure, which traces the square of the gas density in the clus- 
ter. The emission measure peaks on the subcluster core and there is 
clearly a sharp density drop at the core's leading edge to the SE. The 
emission measure declines more gradually through the ram pres- 
sure stripped tail of gas to the NW of the core and there is a plume 
of emission extending to the SW perpendicular to the merger axis 
(Fig. [2] upper left). There is no obvious second peak in the emission 
measure map corresponding to the primary cluster core, although 
there are clumps of more dense material in the approximate posi- 
tion of the main collision site. 

The temperature map (Fig.|4] upper right) shows that the dense 
subcluster core contains the coolest gas in the cluster, down to 
1.45 ±0.07 keV. The steep density drop of the SE edge of the 
cluster core corresponds to a sharp increase in the gas temperature 
from ~ 2 keV to ~ 6 keV. This is a cold front or contact discon- 
tinuity created as the subcluster's dense cool core travels through 
the hotter, diffuse outskirts of the primary cluster ( [Markevitch et al.[ 
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Figure 2. Upper left: Exposure-corrected image in the 0.3-7.0 keV energy band smoothed with a 2D Gaussian (T = 1.5 arcsec (North is up and East is to the 
left). Upper right: Unsharp-masked image created by subtracting images smoothed by 2D Gaussians with a = 5 and 20 arcsec and dividing by the sum of the 
two images. Point sources were removed before unsharp-masking. Note that a point source has been removed at the NE end of the bow shock. Lower left: 
same as upper left but showing the subcluster core (units photons cm^^ s^' pixel"' ). The AGN nucleus is marked by the white dashed circle. Lower light: 
Subaru R-band image of the galaxy distribution in Abell 2146 with the X-ray gas contours overlaid. Concentrations of galaxies are marked with dashed circles. 



|2000[ IVikhlinin, Markevitch & Murray|[2001a[ >. The temperature 
increases more steadily through the ram pressure stripped tail up to 
~ 6 — 7 keV. The plume extending to the SW is cooler than its sur- 
roundings with a temperature of 5 — 6 keV. There does not appear 
to be a similar cool, dense plume structure on the other side of the 
subcluster tail to the NE. 

The hottest gas in the cluster, at temperatures from 12 — 
15 keV, is located at the core collision site to the NW. The NW 
edge of this high temperature region is the upstream shock and cor- 



responds to a drop in the gas density shown by the emission mea- 
sure map. The bow shock is also visible as a peak in the tempera- 
ture map to the SE but this is not as clearly shown by the selected 
binning. 

The projected 'pressure' map (Fig.|4] lower left) was produced 
by multiplying the square root of the emission measure and the tem- 
perature maps. The 'pressure' map reveals the extent of the shock 
heating and compression along the merger axis and the sharp drops 
in pressure to the NW and SE correspond to the upstream and bow 
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Figure 3. Exposure-corrected image of the subcluster core in the 0.3- 
7.0 keV energy band smoothed with a 2D Gaussian fj = 1 arcsec (North 
is up and East is to the left; units photons cm^^ s^' pixel^' ). The blue 
contours were produced from the HST F606W archival image of the galax- 
ies (blue), including the BCG, and the white cross shows the location of 
the VLA 1 .4 GHz point source marking the nucleus (NRAOA'LA Archive 
Survey). 

shock fronts, respectively. The SW plume appears to be in pressure 
equihbrium with its surrounding environment and this is explored 
further in section [2^.3.2| Although there are regions of high pressure 
that could be associated with the core collision site, it was difficult 
to determine this location exactly and the formation of the plume 
indicates it could instead be closer to the subcluster tail (section 

The signal-to-noise ratio was increased to 70 (~ 5000 counts) 
to create the metallicity map. The errors on the metallicity values 
are approximately ±0.08 Zp, but increase to 0.1 — O.lSZrr^ for the 
highest temperature shock-heated regions where the ICM is almost 
completely ionised and there is minimal line emission. The metal- 
licity peaks in the subcluster core, ahead of the BCG, at ~ 0.9 Ztj 
and is approximately constant elsewhere at --^ O.4Z0. The appar- 
ent drop in the metallicity inside the subcluster core is an artifact 
caused by the use of a single temperature spectral model where 
the cluster gas has multiple temperature components ( [Buote &] 
|Canizares|[l994[ |Buote||2000l l. If we fit the spectra from this re- 
gion with a two temperature model, the best-fit metallicity value is 
O.SS^QQgZQ, which is consistent with the average for the subclus- 
ter core. Two temperature model fits are discussed in more detail in 
section|4] 

The metallicity drops sharply across the SE edge of the sub- 
cluster core showing clearly the difference in origin of the gas on 
either side of this contact discontinuity. It is not clear why the post- 
shock gas ahead of the subcluster core has a lower average metal- 
licity. By fitting the spectra from the three low metallicity regions 
ahead of the subcluster core together, we find a best-fit value of 
only O.15±O.O6Z0 compared to the ambient value of --^ O.4Z0. 
Although the gas temperature is high here and the metallicity is 
more difficult to constrain, we note that the metallicity for the 
shock-heated gas behind the upstream shock has a metallicity con- 



sistent with the average at 0.5 ±O.2Z0. There was also no evidence 
from spectral fitting for a second temperature component in this 
region. We consider the possible impact of non-equilibrium ion- 
ization behind the shock fronts in section [34] The metallicity also 
drops rapidly behind the subcluster core from 0.9 ± 0.08 keV to 
O.S^q q^Zq in the region immediately behind the AGN. There is 
no evidence for a metallicity gradient in the ram pressure stripped 
tail suggesting that the metal enriched material is pulled off the sub- 
cluster core in clumps which don't efficiently mix with the ambient 
ICM. 

2.3.2 Radial profiles 

Several key sectors of the merging cluster are identified in Fig. |5] 
and used to investigate these structures in greater detail. We ex- 
tracted projected surface brightness, temperature and metallicity 
profiles in these regions using a minimum of 50 source counts per 
radial bin for the surface brightness and 2000 counts per bin for the 
spectral fitting. The temperature and metallicity were determined 
by fitting a single temperature model to each extracted spectrum as 
in section [2^.3.1| Fig.|6]shows the radial profiles for the SE and NW 
sectors positioned along the merger axis. The SE sector (Fig.jSj is 
centred on the subcluster core and includes the leading edge of the 
core and the bow shock. The NW sector was selected to cover the 
subcluster tail, the collision site and the upstream shock. Note that 
the centre of each sector was not the centre of curvature of each 
shock and so these features appear smoothed here compared to the 
later analysis in section [3TT| 

The surface brightness edges corresponding to the shock 
fronts and the contact discontinuity ahead of the subcluster core 
are clearly visible in Fig.|6]but we also note two additional edges 
in the NW sector immediately behind the subcluster core at 26 kpc 
and at the end of the subcluster tail at 240 kpc. The surface bright- 
ness drop immediately behind the subcluster core is also accom- 
panied by an increase in the gas temperature from 3.2 ±0.1 keV 
to 5.0 ±0.3 keV. Interestingly, the AGN and centre of the BCG 
is located almost exactly at the position of this surface brightness 
and temperature edge. The AGN was excluded from the spectral 
fits using a region of radius 2" and therefore did not contribute to 
the temperature increase at this location. It is more likely that this 
contact discontinuity has been created by the convergence of the 
surrounding ICM behind the dense subcluster core. The metallicity 
immediately behind the core appears consistent with the core value, 
although with large errors. The metallicity then decreases to O.4Z0 
and is approximately constant through the subcluster tail. Metal- 
rich material that is ram pressure stripped from the subcluster core 
may build up in the region immediately behind the core where the 
flow of the ambient ICM converges. 

The surface brightness edge at the end of the subcluster tail 
(at 240 kpc) is broad, covering ~ 80 kpc in radius, and has a rel- 
atively shallow surface brightness gradient that does not appear to 
be characteristic of a contact discontinuity. The temperature gradi- 
ent increases from 6.7 ± 0.2 keV to 9. I^q g keV across this region. 
However, there is likely to be a significant component of projected 
emission on this region which will increase the apparent temper- 
ature. Deprojection routines usually assume spherical symmetry 
which is not a good assumption for this region of Abell 2146 where 
the cluster emission is strongly elongated along the merger axis and 
the surface brightness gradient is relatively shallow. This surface 
brightness edge marks the region where 6 — 7 keV material from 
the end of the subcluster tail could be mixing with gas from the 
primary cluster. 
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Figure 4. Upper left: projected emission measure per unit area map (units are logio cm^ arcsec^ ). The emission measure is the XSPEC normalization of 
the MEKAL spectrum K = EI / {Ax 10''*;rD^(l + z)^), where EI is the emission integral EI = J ncrtfjdV. The approximate position of the core collision site 
is labelled. Upper right: projected temperature map (keV). Note that the white circle behind the cool core is the excluded central AGN. Lower left: projected 
pseudo-pressure map (keVcm^^^- arcsec"-). Lower right: projected metallicity map (Z,; ) generated using larger spectral bins with ~ 5000 counts. The 
excluded point sources are visible as small white circles. North is up and East is to the left. 



Fig. |7] shows the surface brightness and temperature of the 
gas in the SW plume and a comparison no-plume sector on the 
other side of the subcluster tail (see Fig.|5j. These two sectors were 
also each compared with two neighbouring sectors positioned ei- 
ther side to the NW and SE. Fig. |7] (left) shows the surface bright- 
ness enhancement of the plume over the surrounding sectors from 
~ 50 — 250 kpc. The electron pressure is possibly higher in the sec- 
tor to the SE of the plume as the temperature is higher here but 
the errors are large, otherwise the plume is in pressure equilibrium 
with the surroundings. There is no comparable surface brightness 



increase in the no-plume sector. The increase in temperature in the 
NW no-plume sector is probably due to the inclusion of a section 
of the upstream shock in the region analysed. The metallicity is 
approximately constant at OAZq, similar to the ambient ICM, for 
both plume and no-plume sectors. 

We suggest that the plume is likely to be the remnant of the 
primary cluster core which has been pushed forward and laterally 
by the impact of the subcluster core. There does not appear to be a 
symmetric structure extending from the merger axis to the NE but 
this can be explained if the two clusters collided with a non-zero 
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Figure 5. Image showing the sectors used to produce the radial profiles in 
Fig.|6]and Fig.|7]with the surface brightness edges marked by dashed white 
lines. 



impact parameter (see eg. Roettiger, Stone & Mushotzky||1998l 



IRicker & Sarazin|2001[ [Poole et al.|2006^ . Simulations of off-axis 
mergers indicate that if the subcluster passed to the North of the 
primary cluster centre, at a distance of order the core size, material 
from the primary core would be ejected perpendicular to the merger 
axis in the SW direction, as observed in Abell 2146. A strongly off- 
axis collision, comparable to the lateral extent of the shock fronts, 
seems unlikely however as the two shock fronts are both close to 
symmetric about the merger axis and the primary cluster core has 
been strongly disrupted by the merger. Detailed simulations of the 
merger in Abell 2146 will be needed to improve these rough limits. 



3 THE SHOCK FRONTS 
3.1 Surface brightness profiles 

Fig.[8]shows the sectors selected for the analysis of the bow and up- 
stream shock fronts. The sectors were positioned according to the 
centre of curvature of the edge and cover an angular range where 
the shock front is well-defined. These sectors were divided into ra- 
dial bins initially of 1 arcsec width but this was increased to 1.5, 
2.5 and 5 arcsec at larger radii to ensure a minimum of 50 source 
counts in each bin. Point sources were excluded from the analy- 
sis and the background was subtracted using the normalized blank 
sky background datasets detailed in section [2T| To maximise the 
signal-to-noise ratio, the surface brightness profiles were restricted 
to the energy range 0.5 — 4.0 keV and we used only the Chandra 
observations taken with ACIS-1, which have a lower background 
than ACIS-S. Fig. |9] shows the final background-subtracted X-ray 
surface brightness profiles for the two shock fronts. We fit the sur- 
face brightness edge associated with each shock with a model for 
a projected spherical density discontinuity, assuming that the edge 
is symmetric about the merger axis (eg. [Markevitch et al.||2000| 
|2002[[Owers et al.|2009| l. The model radial gas density profile con- 
sists of a power law on either side of an abrupt density jump where 
the free parameters are the slopes and normalizations of the power 
laws and the radial position of the density jump. The corresponding 
emission measure profile was projected onto the sky and fitted to 
the observed surface brightness profile over a radial range selected 



to exclude unrelated core structures but extending to large radii for 
correct projection. 

Fig. |9] shows that this model provides a reasonable fit to both 
the bow and upstream shock fronts. The amplitude of the density 
jump at each shock was calculated from the square root of the ra- 
tio of the power law normalizations (eg. |Owers et al.|2009] l with 
a small correction for the observed temperature difference across 
the edge. Following [Landau & Lifshitz| ( [1959^ , we applied the 
Rankine-Hugoniot shock jump conditions across each shock front 
to calculate the respective Mach numbers, M = v/cj, where v is 
the velocity of the preshock gas with respect to the shock surface 
and Cs is the velocity of sound in that gas. The Mach number was 
calculated from the density jump. 



M = 



y+1- ''Vp,(r-i) 



1/2 



(1) 



where pi and P2 denote the gas density upstream and downstream 
of the shock respectivel}]^ We assume the adiabatic index for a 
monatomic gas, 7=5/3. 

At the bow shock, the density drops by a factor Pi/pi = 2.5 ± 
0.2 which gives a Mach number M = 2.3 ± 0.2. For the upstream 
shock, the density decreases by P2/P1 = 1.9 ±0.2 producing M = 
1.6±0.1. 

Even in the raw counts image shown in Fig.[T|it is clear that 
the upstream shock in Abell 2146 has a very narrow edge separating 
the pre- and postshock gas. We estimated the width of the shocks 
by smoothing the best-fit density discontinuity model with Gaus- 
sian functions of varying widths, CT, and fitting the smoothed mod- 
els to the surface brightness profile across the shock front. For the 
upstream shock, the best-fit model with zero width has = 98.4 
for 85 degrees of freedom. This is reduced to = 88.9 for 84 de- 
grees of freedom for a smoothed model with a width CT = 6'^^ ^pc 
(95% errors). Fig. |9] shows that the bow shock is broader than the 
upstream shock with a best-fit width of 12^^ kpc (95% errors). 

The shock widths can be compared with the Coulomb mean 
free path of the electrons and protons on both sides of the shock 
front. Following [Vikhlinin, Markevitch & Murray|j2001a| l, we es- 
timate the mean free path for particles crossing from the postshock 
gas into the preshock gas, Ain^om^ which is the main source of 
diffusion across the edge. For the bow shock, the mean free path 
■^n^out = 21 ± 3 kpc and for the upstream shock it is 23 ± 5 kpc. If 
Coulomb diffusion is not suppressed across the shock edge, we ex- 
pect the shock to have a width of several times the mean free path. 
The upstream shock is significantly narrower than the estimated 
mean free path, which suggests there is a significant suppression of 
transport processes across this edge consistent with a collisionless 
shock. 

The bow shock appears broader however the measured width 
of the shock fronts can be affected by any deformations in the front 
shape, which smear the edge in projection. It is therefore difficult 
to determine if the bow shock front is intrinsically wider. The esti- 
mated width of the Bullet cluster bow shock is ~ 35 kpc, although 
|Markevitch,(2006) find this is only marginally preferred over a zero 
width shock. It is however surprising that the upstream shock front 



^ Note that eq. 2 in [Russell et al.[|2010| , relating the Mach number to the 
temperature jump, is applicable only for strong shocks. Given the large eiTor 
on the temperature values in jRussell et al.[j2010) this had a negligable effect 
but as the temperature errors are much smaller for this deeper Chandra 
observation a similar assumption is not made here. 
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Figure 6. Surface brightness, temperature and metallicity profiles for the SE (left) and NW (right) sectors shown in Fig. [3] The surface biightness profiles 
cover the energy range 0.5 — 7 keV and edges, corresponding to cold fronts or shocks, are marked by dashed lines. Note that the sectors were not selected to 
correspond to the shock edges (see Fig. |l 1) . 



should be the narrowest when it is propagating through the remain- 
ing infall from the subcluster's passage. 

Fig. [To] shows the variation in the surface brightness profile 
of each shock front with position angle. The bow shock appears to 
be slightly stronger in the central sector but the lower number of 
counts in this sector, particularly in the preshock region, increases 
the error on the density jump value and this is therefore not a sig- 
nificant result. Both the bow and upstream shock fronts appear to 
be consistently narrow across their respective lengths of ~ 500 kpc 
and ~ 440 kpc. 



3.2 Spectral analysis of the shock fronts 

The temperature jump across the shock front can be used to de- 
rive an independent measure of the Mach number and to calcu- 
late the shock velocity to determine the time since core passage of 
the subcluster Fig. [TT] shows the projected temperature, projected 
metallicity and deprojected electron density profiles for the bow 
shock and upstream shock sectors (see Fig.[8](. These radial profiles 
have the same radius of curvature as each of the two shock fronts 
and therefore clearly show the temperature and density decreases 
across the shock edges. We calculated the deprojected temperature 



profile across each shock front using PROJCT in XSPEC. Deprojec- 
tion routines, such as PROJCT, generally assume that the cluster is 
spherically symmetric, which is a reasonable assumption for a re- 
laxed cluster but can be problematic for major mergers. However, 
for the sectors across the shock fronts, the cluster appears approxi- 
mately circular on the sky and the sharp increase in surface bright- 
ness across the edge reduces the impact of projected emission. The 
large size of the radial bins used to measure the temperature values 
is likely to be a more significant source of error in determining the 
temperature jump across each shock. 

We find that the deprojected temperature decreases by a fac- 
tor of T2/T1 = 1.8 ± 0.3 at the bow shock and T2/T1 = 2. l+'J-^* at 
the upstream shock. The Rankine-Hugoniot shock jump conditions 
directly relate the temperature jump to the density jump. We can 
therefore calculate the expected density jump and Mach number of 
each shock (eq. [TJ from the observed temperature jump as an in- 
dependent verification ( [Landau & Lifshitzl|1959| l. We calculate a 
Mach number M = l.S+Jj j for the bow shock and M = 2.0 ± 0.3 
for the upstream shock. The Mach number for both shock fronts is 
therefore consistent within the errors for the calculations using both 
the temperature and density jumps. The temperature jump provides 
a less accurate measure of the Mach number as the errors are greater 
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Figure 7. Surface brightness and temperature profiles for the plume and no-plume sectors shown in Fig. [5] The NW and SE profiles refer to additional, 
comparison sectors positioned either side of the plume and no-plume sectors. 
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Figure 8. Unsharp-masked image (as in Fig. [2] upper right) showing the 
sectors used to produce surface brightness profiles for the shock fronts. The 
sectors were centred on the centre of curvature for each shock front. 



than for the density and larger radial bins must be used to measure 
the temperature so we cannot resolve the jump across the shock 
accurately. 

By using the bow shock speed and estimating the distance be- 
tween the subcluster core and the primary cluster centre, we cal- 



culate the time since the subcluster passed through the primary 
cluster core. The shock speed, v, is calculated by multiplying the 
shock Mach number and the sound speed in the preshock gas 
Cj = (YkBTi/niiiil)^/'^, where Ti is the preshock gas temperature 
and fl = 0.6 is the mean molecular weight in the ICM. For the 
bow shock in Abell 2146, we use the Mach number calculated 
from the density jump M = 2.3 ± 0.2 and calculate the preshock 
sound speed to be Cj = HVOI'^'q" km s^' , giving a shock velocity 
V = 2700^3QQ km s^' . However, even in the new, deep Chandra 
observations, it is not possible to conclusively identify an X-ray 
peak associated with the primary cluster core (Fig.|4] upper left). 
Simulations suggest that if the primary cluster halo has a low con- 
centration, the core will not survive a major collision and the X-ray 
peak will be destroyed (eg. |Mastropietro & Burkert|2008p . We have 
therefore assumed that the collision between the two cluster cores 
occurred at a peak in the pressure map between the subcluster tail 
end and the upstream shock (Fig. |4] left). The distance between 
the subcluster core and the estimated location of the core collision 
site is then --^ 350 kpc and the time since core passage is therefore 
--^ 0.1 — 0.2 Gyr. This is only a rough estimate of the timescale as 
the subcluster velocity is likely to be significantly lower than the 
shock velocity ( [Spring el & Farrar||2007|l and the location of the 
collision site may be closer to the subcluster tail and the plume. 
However, it is clear from the detection of the shock fronts and the 
SW plume structure that this merger is observed soon after core 
collision. 

Fig. [TT] (left) also shows that there is a sharp change in the 
gas properties across the cold front at the leading edge of the sub- 
cluster cool core. The temperature increases by a factor of 3.7+°-^ 
(discussed further in section [4~2^ whilst the metallicity drops from 
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Figure 9. Left: Background-subtracted surface brightness profile across the bow shock in the energy band 0.5-4.0 keV overlaid with the best fit projected 
density discontinuity model with no smoothing (solid lines) and with 20 kpc smoothing, corresponding to the mean free path X^^^^m (dashed lines). Right: 
same as left but for the upstream shock. 




Figure 10. Background-subtracted surface brightness profiles of the bow shock (left) and upstream shock (right) each divided into three sectors showing the 
variation in the shock width and strength with position angle. The best-fit projected discontinuity models generated for the full width of the shock fronts in 
section lxTl are shown overlaid. 



0.89lo Q5 ^0 inside the subcluster core to O.SO+Q og Zq in the post- 
shock gas. The sharp metallicity drop clearly demonstrates the dif- 
ferent origins of the ICM on either side of the contact discontinuity. 
The high metallicity core from the subcluster is travelling through 
the ICM in the outskirts of the primary cluster, which has a much 
lower metallicity. In comparison, the upstream shock sector (Fig. 
[TT|right) has an approximately constant metallicity of 0.3 — O.4Z0, 
as it consists predominantly of primary cluster gas. There may be 
more metal-enriched ram pressure stripped material from the sub- 
cluster in this sector, however the material that was stripped from 



the core earlier is likely to have a lower metallicity than the current 
peak and the subcluster tail shows no evidence of a strong metal- 
licity gradient (Fig.|4]l. 

3.3 The establishment of electron-ion equilibrium 

Cluster merger shock fronts present a unique opportunity to inves- 
tigate the electron-ion equilibration time in the magnetized ICM 
by mapping the postshock electron temperature. The Rankine- 
Hugoniot shock jump conditions can be used to calculate the post- 
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Figure 11. Radial profiles through the bow (left) and upstream (right) shock sectors (shown in Fig. [sj showing the projected electron temperature (upper), 
projected metallicity (centre) and deprojected electron density (lower). There are two sets of points showing a finer (filled points) and a broader (open points) 
radial binning. The electron density model (solid line) shown corresponds to the best-fit to the surface brightness profile across each shock (section fxT) . The 
vertical long-dashed line shows the location of each shock front and the vertical short-dashed line shows the position of the cold front on the leading edge of 
the subcluster core. 



shock temperature for the ICM electrons and ions once they reach 
equilibrium after the passage of the shock (eg. [Landau & Lifshitz] 
|1959j ). However, the fraction of a shock's kinetic energy that is ini- 
tially transferred to the thermal and cosmic-ray populations of the 
electrons and ions remains an open question. 

A shock front propagating through a coUisional plasma heats 
the ions dissipatively in a shock layer that has a width of order the 
coUisional mean free path. The electrons have a much higher ther- 
mal velocity and are not strongly heated by merger shocks. They 
are compressed adiabatically and subsequently equilibrate with the 
ions according to the Coulomb coUisional timescale l ,Spitzer,1962) 
given by. 



layer several orders of magnitude thinner than the mean free path 



j{e,p)^6.2x 10** yr 



T \3/2 
108 K 



10- 



(2) 



where Tg is the electron temperature and is the electron density 
(seeeg. [Sarazin|1988) . 

However, shocks in a magnetised plasma, such as the ICM, 
are likely to be coUisionless. Observations of the solar wind shocks 
found that the electron and ion temperature jump occurs in a shock 



(eg. Ness, Scearce & Seek U964[ [Montgomery, Asbridge & Bame| 
[1970[ [Hull et al. |2001 1. The coupling of particles to electric and 
magnetic fields produces interactions which have dissipation scale 
lengths much shorter than the ordinary collision mean free path. 
The plasma waves producing these interactions and instabilities 
affect ions and electrons differently due to the large difference in 
mass (for a review see eg. Tidman & Krall|1971[ [Friedman et aT] 
[1971[ l. We might therefore expect to find an electron heating rate 
shorter than the Coulomb coUisional timescale behind a cluster 
merger shock. 

Measurements of the postshock temperatures in supernova 
remnants (eg. jR akowski 2005 , Raymo nd & Korreck[ [2005[ 
[Ghavamian, Laming & Ra kowski 2007) and heliospheric shocks 
(eg. [Schwartz et al.[[1988[ [Russell 20(35} show that in most cases 
the electrons are heated less than the protons, Tf/T, < 1, in regions 
close to the shock front. However, these measurements cannot de- 
termine the timescale of subsequent equilibration, teq{e,p), as this 
corresponds to a linear scale of several A.U. for planetary bow 
shocks and is comparable to the age of the remnant for supernova 



© 0000 RAS, MNRAS 000, 000-000 



The merging cluster Abell 2146 13 



17.5 - 



CoUisional equilibration (tfm) 

Instant equilibration (sfm) 



12.5 



I 




800 1000 
Radius (kpc) 

Figure 12. The re-establishment of electron-ion equilibrium behind the bow 
shock. The projected electron temperature profile is overlaid with model 
predictions (with 1<T error bands) for instant equilibration (red) and adi- 
abatic compression followed by collisional equilibration (blue). The open 
and filled data points show narrower and broader binning of the tempera- 
ture profile. 
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Figure 13. The re-establishment of electron-ion equilibrium behind the up- 
stream shock. The projected electron temperature profile is overlaid with 
model predictions (with Iff error bands) for instant equilibration (red) 
and adiabatic compression followed by collisional equilibration (blue). The 
open and filled data points show nan'ower and broader binning of the tem- 
perature profile. 



shocks. For cluster merger shock fronts the equilibration timescale 
corresponds to a distance of hundreds of kpc, so while clusters are 
too distant for us to resolve the shock layer, we can study the sub- 
sequent equilibration of the ICM constituents. 

The bow shock in the Bullet cluster provided the first measure- 
ment of the electron-ion equilibration timescale in the ICM. |Marke-| 
|vitch| j2006| compared the observed electron temperature profile 
across the shock front with two models for equilibration. The in- 
stant equilibration model predicts that the electrons are strongly 
heated at the shock front and the electron temperature rapidly in- 
creases to the postshock temperature, similar to the ion tempera- 
ture. The collisional model predicts an adiabatic compression of 
the electron population at the shock and a subsequent slower equi- 
libration with the ions on a timescale determined by Coulomb col- 
lisions. The observed temperature profile for the Bullet cluster sup- 
ported instant equilibration, indicating that electrons were rapidly 
heated at the shock front on a timescale faster than Coulomb colli- 
sions. However, the postshock temperature in the Bullet cluster is 
very high (~ 20 — 40 keV) compared to the Chandra energy band 
and therefore difficult to constrain. 

Although the postshock temperatures in Abell 2146 are lower 
than the Bullet cluster, the shock Mach numbers are also lower, 
which reduces the separation between the collisional and instant 
equilibration models. Figs.|I2|and|13|show the observed and model 
projected electron temperature profiles behind the bow shock and 
upstream shock, respectively. The observed profiles were generated 
using single temperature model fits to spectra extracted from two 
different sets of radial bins. The narrower radial bins were selected 
to contain a minimum of 2000 counts and the broader bins have a 
minimum of 4000 counts. The postshock radial bin closest to the 
shock boundary was positioned to overlap the fainter preshock gas 
and ensure minimal contamination of the neighbouring preshock 
gas radial bin. 

The models were generated by using the observed preshock 
electron temperature and the density jump in the Rankine-Hugoniot 
shock jump conditions to predict the postshock gas temperature. 
For the instant equilibration model (a single fluid model), the elec- 
tron temperature and the ion temperature jump at the shock front 
to the postshock gas temperature predicted by the shock jump con- 
ditions. This temperature model was then projected along the line 
of sight, using the best-fit surface brightness model determined for 
each shock front in section [3TT| to produce the final projected tem- 
perature model. The projection smoothes the shaip jump in temper- 
ature at the shock front. Although the projection assumes a constant 
preshock temperature, the steep decline in the surface brightness 
across and ahead of each shock front (Fig. |9]l meant the effect of 
any cooler gas beyond the observed temperature points is negliga- 
ble. 

For the collisional equilibration model (a two fluid model), the 
electron temperature increases at the shock front according to the 
adiabatic compression of the particles. 



(3) 



Coulomb collisions then subsequently equilibrate the electron and 
ion temperatures at a rate 



dr^ 



(4) 
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We require that the total kinetic energy density is conserved so that 
the local mean gas temperature is given by 

_ njTj + neTe _ 7] + 1.17; 
ngas ' 2.1 • 

The mean gas temperature is constant and can be calculated from 
the Rankine-Hugoniot shock jump conditions. Therefore, although 
X-ray observations do not currently have the spectral resolution 
required to measure the ion temperature directly, we can use this 
requirement to determine the immediate postshock ion tempera- 
ture. The postshock ion temperature is significantly higher than 
the adiabatically-compressed electron temperature and for the col- 
lisional model we assume that the electrons equilibrate with the 
hotter ions according to the Coulomb coUisional timescale, teq. Us- 
ing eq. |2] we calculate this timescale for the postshock electrons. 
By multiplying teg by the shock velocity in the postshock gas (bow 
shock Vps = 1100± 100 km s^', upstream shock Vps = 1300 ± 
200kms-i), 

we determine the distance behind the shock where 
equilibration is reached. Electrons at this position were heated by 
the shock 1,,^ years ago. By integrating eq.|4] and using eqs.|2]and|5] 
the electron temperature as a function of distance behind the shock 
was determined analytically (see eg. |Fox & Loeb|199'7l [Ettori &| 
|Fabian|1998[ l. This collisional model was then also projected along 
the line of sight to produce the final model for comparison with the 
observed projected electron temperature. 

The instant and collisional models were projected by deter- 
mining the model electron temperature in small volumes, dV, along 
the line of sight for a particular annulus and calculating the corre- 
sponding emission measure using the best-fit density discontinuity 
model (section [3?7| (. For each annulus, the emission measures for 
dV with similar temperatures were summed together using a set 
temperature binning with fine resolution (~ 0. 1 keV). To determine 
the projected temperature in this annulus, we produced fake spectra 
in XSPEC using multi-component MEKAL models with the temper- 
ature of each component set to the midpoint of each temperature 
bin and the normalization set to the corresponding summed emis- 
sion measure. The metallicity was fixed to 0.4 Zq and responses 
for the appropriate detector region were used. The MEKAL model 
components were also combined with a PHABS absorption compo- 
nent set to the Galactic column density. These fake spectra were 
then fitted with single absorbed MEKAL models, with fixed metal- 
licity and column density, to determine the final model projected 
temperature. 

The main sources of error for these models were the measure- 
ment of the preshock temperature and, to a lesser extent, the density 
jump and preshock electron density. We used a Monte Carlo tech- 
nique to determine the uncertainties in the model projected tem- 
perature profiles. We repeated the model calculation and projection 
1000 times, each time using new values of the preshock temper- 
ature, density jump and preshock density based on Gaussian dis- 
tributions. The output models are the median of this process and 
the IcT errors are calculated from the 15.85 and 84.15 percentile 
spectra. 

Fig.[T2]shows the observed and model projected electron tem- 
perature profiles across the bow shock. The observed temperature 
profile is cut off ~ 150 kpc behind the bow shock front where 
the temperature drops in the subcluster core. Unfortunately, ow- 
ing to the large uncertainty on the preshock count rate in the earlier 
45 ks observation, there were fewer preshock counts in the new 
400 ks observation than anticipated. The lower number of counts 
produced a greater than expected error on the crucial preshock 



temperature, which is the main source of uncertainty for the mod- 
els, and therefore it is not possible to conclusively exclude either 
model. However, the observed postshock electron temperatures for 
the bow shock appear lower than predicted by the instant equili- 
bration model and favour the collisional equilibration model. We 
therefore conclude that collisional equilibration cannot be ruled out 
for cluster merger shocks. 

The preshock temperature could potentially be higher than the 
conservative estimate of 5.3i'iQ ^ keV used in this analysis. The tem- 
perature bin closest to the bow shock front has a slightly higher 
temperature of 6.4^Qg keV. This increase in temperature could be 
due to contamination from the higher temperature and density post- 
shock region. However, the neighbouring postshock temperature 
bin was positioned to significantly overlap the preshock region by 
20 kpc (6 arcsec) and minimise this effect. It is therefore more likely 
that there is a gradient in the preshock gas temperature which would 
increase the postshock temperature of both models and strengthen 
our conclusion that collisional equilibration is possible for cluster 
merger shocks. 

Measurement of the equilibration timescale was more difficult 
for the upstream shock as the Mach number was lower, reducing the 
separation between the instant and collisional models. In addition, 
the postshock electron temperatures are ~ 5 keV higher compared 
to the bow shock, which increased the uncertainty on the values. 
Fig.[T3]shows that the electron temperature increases rapidly in the 
postshock region to 15l[ keV and then drops to ^ 10 keV approx- 
imately 100 kpc behind the shock front. Although the errors are 
larger, the electrons behind the upstream shock appear to equili- 
brate faster than those at the bow shock. However, the upstream 
shock is also propagating through the primary cluster outskirts, 
which have been disturbed by the subcluster's passage. There is 
likely to be significant substructure and additional shock heating 
in the region of the upstream shock front. Ram pressure stripped 
material from the subcluster tail is likely to be the cause of the 
rapid decline in temperature ~ 70 kpc behind the shock. We there- 
fore place less weight on the conclusions on equilibration behind 
the upstream shock as the situation here is more complex and the 
errors on the electron temperature are greater. 

3.4 Non-equilibrium ionization 

In the previous section, we have shown that the postshock ICM 
significantly deviates from thermal equipartition between electrons 
and ions but our analysis has so far assumed that the ICM is in 
ionization equilibrium. However, while the temperature of the ICM 
increases rapidly in the shock layer, the ionization state of the ions 
still reflects the preshock temperature and the ICM will be under- 
ionized compared to the equilibrium case. The ionization balance 
will be recovered by collisions on a timescale of --^ 10^ yr for an 
electron density --^ 10^^ cm^-^ but until this is achieved there will 
be more ionizations than recombinations. Simulations of the outer 
regions of clusters and of merging systems have considered the ef- 
fects of non-equilibrium ionization, in particular the effect on the 
intensity ratios of the Fe K a lines (eg.|Yoshikawa & Sasa ki|2006[ 
[Akahori & Yoshikawa 2010 201 1 , Wong, Sarazin & Ji 20flJ. We 
have therefore considered whether non-equilibrium ionization will 
produce an observable effect in the postshock ICM of Abell 2146. 

For a temperature of '-^ 9 keV (10** K) behind the bow shock, 
the typical density-weighted timescale for the ions to approach 
ionization equilibrium is T --^ 4 x lO'^ cm^-' s for Fe and signif- 
icantly shorter for other important ions in the ICM (eg. [Smith &] 
|Hughes|2010| l. Using the postshock electron density, 1 .99 ± 0.09 x 
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10~ cm~^ , and the velocity of the shock in the postshock gas (sec- 
tion |3.3[ l, we estimate that the Fe ions reach equilibrium ionization 
^ 50 kpc behind the bow shock. Therefore, only the first narrow 
radial bin behind the bow shock (and also for the upstream shock) 
could be affected by non-equilibrium effects. 

We fit each of these radial bins, immediately behind each 
shock front, with an absorbed non-equilibrium ionization NEI spec- 
tral model in XSPEC (eg. [Hamilton, Sarazin & Chevalier||1983| 
|Borkowski, Lyerly & Reynolds 200 l l. The NEI model indicated 
a density-weighted ionization timescale of T > 4 x lO'-^ cm^^ s 
behind each shock front suggesting that there was no measurable 
signature of non-equilibrium ionization in Abell 2146. This is a 
difficult measurement given the low spectral resolution at high en- 
ergies, covering the Fe XXV K a and Fe XXVI K a lines ( |Akahori| 
|& Yoshikawa||20l"0l|2011[ l, and the low photon count rates at the 
shock front, necessitating large radial bins relative to the ionization 
equilibration timescale. 



4 DISRUPTION OF THE SUBCLUSTER CORE 

Fig. |14| shows the detailed structure in the subcluster core and ram 
pressure stripped tail. The front edge of the subcluster core ap- 
pears roughly spherical but the sides are strongly sheared by the 
surrounding shocked gas. The southern edge of the subcluster core 
appears to be narrow with no sign of the disruption and stripping 
which affects the eastern and northern edges (eg. Fig.|2]l. The tem- 
perature gradients on the SW and NE sides of the subcluster core 
are very different. There is clearly a more graduated increase in 
temperature along the NE edge where the core is breaking up. 

Fig.[T5]shows that the filament of gas running along the south- 
em edge of the core, and apparently ending at the AGN, contains 
the coldest gas in the cluster. Using a single temperature spectral 
model (section|23l, we determined that the X-ray gas temperature 
drops to 1 .48ly gg keV in this filament. This cool X-ray gas fil- 
ament is also detected as a narrow, coherent filament in Ha and 
[N II] observations ( [Canning et al.|20TT} . We fitted the spectrum 
from this region (shown in Fig. 1 15} with a multi-temperature model 
to determine if there is a significant amount of multi-phase gas de- 
tectable in the X-ray. Following |Sanders et al.| (|2004'), we used mul- 
tiple absorbed MEKAL components with temperatures fixed at 0.5, 
1, 2 and 4 keV, common metallicities and free normalizations. We 
tested an additional 8 keV component but found that the best-fit 
normalization was consistent with zero. The majority of the gas in 
the filament is at around 2 keV but we detect significant amounts of 
cooler gas. The fraction of the emission measure of gas at 1 keV, 
with respect to the total emission measure, is 15 — 20% and at 
0.5 keV is 5 — 10%. Both of these lower temperature components 
are detected at above 3(7. The best-fit metallicity for this multi- 
temperature model is 0.8±0. IZg. It is therefore likely that the 
low metallicity observed in the subcluster core in Figs.[4]an d[6|is a 
bias caused by the use of only a single temperature model ^Buote| 
|& Camzares[T994l|Buote|2000) . 

Using an absorbed MEKAL + MKCFLOW model, with the 
lower temperature of the MKCFLOW component ( [Mushotzky &| 
|Szymkowiak|1988| > fixed to 0. 1 keV, the upper temperature com- 
ponent tied to the temperature of the MEKAL model and the metal- 
licities tied together, we determine an upper limit on the mass depo- 
sition rate for the filament region of 40 Mq yr^' , which is cooling 
out of the X-ray band (see eg. |McNamara et ar[2 006 , Raff erty et al.| 
|2006| l. The upper limit on the mass deposition rate for the whole 
of the subcluster core is 50Mq yr^^ , therefore the vast majority 



of the cooling in the subcluster core is likely occurring in this fila- 
ment. This is significantly less than the star formation rate (SFR) of 
192M0 yr^' determined from Spitzer observations of the BCG by 
[O'Dea et al. (2()08 l. However the IR emission was not corrected for 
an AGN contribution, which is likely to dominate ( [Canning et al.| 
[201 1| |, and therefore this SFR is an upper limit and the true rate is 
likely to be much lower. 



4.1 Stripping of the NE edge 

Fig. [TS] shows the stripping of cool material from the eastern edge 
of the subcluster cool core. Cool gas is being pulled off the eastern 
edge of the core into a 12" (45 kpc) filament which has a steadily 
increasing temperature from 1.9±0.1 keV to S.l^Jj^ keV along its 
length, although this could be affected by hotter material seen in 
projection. There are several cooler blobs of gas to the North and 
NW of the AGN suggesting this material then breaks off and may 
be directed West into the subcluster' s wake. Warmer 3 — 4 keV gas 
appears to be filling underneath the cool, stripped filament. The 
coolest gas in the subcluster core is towards the leading southern 
edge of the core. 

Simulations of ram pressure stripping of a cool core during 
a cluster merger show strong shearing of the contact discontinuity 
and the formation of a strong vortex inside the core, which trans- 
ports material inside the core to the leading edge and then back 
along the surface travelling with the surrounding flow (eg.! Murr"ay| 
et al. 1993, Bals ara, Livio & 0'Dea| 1994; Heinz et al. 2003). The 
cooler blobs of gas stripped off the eastern edge trace the flow of 
material around the subcluster core and show the flow converges 
behind the core, at the approximate position of the AGN (see sec- 
tion|23]2j. 

In section [2.3.2[ we suggest the subcluster is likely to have 
passed to the North of the primary cluster core, colliding with a 
small but non-zero impact parameter. This trajectory will produce 
a more curved orbit for the subcluster core, trending from a SE 
direction to southern, compared to the case of the head on collision. 
The southern edge is now the leading edge of the subcluster core 
and there is likely to be a greater shearing effect on the eastern edge 
due to the curved trajectory. 

If ram pressure stripping is the dominant mechanism disrupt- 
ing the eastern edge, we expect most of the cool core gas will be 
removed at a radius which satisfies the condition 

P.sV?el>/'cW, (6) 

where Ps is the density of the surrounding ambient ICM, Vfei is the 
relative velocity of the ambient ICM and the subcluster cool core, 
and Pc is the pressure profile of the subcluster core (eg. [Fabian &| 
[Daines[|l991[ [Gomez et al. 2002 1. The non-zero impact parame- 
ter of the merger makes it difficult to determine the velocity of the 
ambient ICM around the eastern edge of the subcluster. We there- 
fore estimate the maximum velocity to be the postshock gas ve- 
locity of 1 100 ± 100 km s" ' . For ;v = 0.0064 ± 0.0002 cm"^ and 
Pq = 0.074 ± 0.006 keV cm^-' , we find that the condition for ram 
pressure stripping is only approximately met along the eastern edge 
of the core although the pressure does not steeply increase here as 
observed inside the southern edge (Fig. [16] lower right). Therefore 
it's likely that the removal of material along this edge is also facil- 
itated by developing Kelvin-Helmholtz instabilities ( [Nulsen|1982[ 
[Inogamov[1999[ l. However, without a more effective estimate of the 
velocity of the ambient ICM around the subcluster core it is difficult 
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Figure 14. Left: projected pseudo-pressure map of the subcluster tail with S/N > 32 in each bin (units keV cm^^''- arcsec"- ). Right: projected temperature 
map (units keV) of the same region with S/A' > 32. The excluded central AGN is visible as a small white circle. The circles mark a region where the temperature 
gradient through the subcluster tail reverses from 6 — 6.5 keV (dashed) to 5 keV gas (solid). 
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Figure 15. Left: projected emission measure per unit area map of the subcluster core with S/N > 15 in each bin (units cm^^ arcsec^^ , see Fig.|4|. Right: 
projected temperature map (units keV) of the same region with S/N > 22. The white dashed region marks the cold gas filament. The excluded central AGN is 
visible as a small white filled circle in both images. 

to calculate the timescale of the developing instabilities and deter- 4.2 Width of the subcluster southern edge 
mine whether the subcluster core will be subsequently destroyed 

(eg. l ViMilinin, Markevitch & Murray '|2001bl| Heinz et £ 112003^ ^" comparison, the SW edge of the subcluster core and tail appears 

to be narrow and stable over a distance of 150 kpc (eg. Fig.Rb. Fig. 

\T6\ (upper right) shows that the surface brightness drops by close 
to an order of magnitude in a distance of only ~ 15 kpc across the 
narrowest point of the southern edge. The eastern edge of the sub- 
cluster core shows a much more gradual decline in surface bright- 



© 0000 RAS, MNRAS 000, 000^000 



The merging cluster Abell 2146 17 




12.5 



10 



7.5 



2.5 



-\- subcluster S edge 
^ subcluster E edge 



20 



50 100 200 
Radius (kpc) 



500 




subcluster S edge 
subcluster E edge 



20 50 100 200 

Radius (kpc) 



500 




0.01 



subcluster S edge 
-\- subcluster E edge 



10" 



20 



50 



100 200 
Radius (kpc) 



500 



Figure 16. Surface brightness (upper right), projected temperature (lower left) and electron pressure profiles (lower right) for two sectors across the southern 
(filled circles, 270 — 315°) and eastern (open circles, 150 — 205°) edges of the subcluster core. An image showing the sectors used is also included (upper left). 
The electron pressure profiles were produced by multiplying the deprojected electron density and the deprojected electron temperature profiles. 



ness with no obvious sharp edge corresponding to a density jump. 
The deprojected temperature, determined using PROJCT in XSPEC, 
increases rapidly across the narrow southern edge of the subcluster 
core by a factor of 3.4+[|'^. The eas in the subcluster core is multi- 
phase and the temperature jump for the cooler 0.5 keV component 
could be higher by a factor of a few. Unless suppressed by at least 
an order of magnitude, thermal conductivity in unmagnetized clus- 
ter gas should rapidl y evaporat e the cool subcluster on a timescale 
of only ~ 10^ yr (eg. Ettori & Fabian 2000, Markevitch et al.,2003, 



|Asai, Fukuda & Matsumoto|2004P 

To determine the density jump across the subcluster south- 
em edge, we fit the surface brightness profile with the model for a 
projected spherical density discontinuity discussed in section [3!T] 
A surface brightness model for the outer cluster emission, a j8- 
model ( [Cavaliere & Fusco-Femiano||1976| > plus a gaussian com- 



ponent to account for the blurred bow shock edge, was added to 
the projected discontinuity model to ensure the projected emission 
was accounted for. Fig. ^] shows the best-fit model for the sur- 
face brightness across the southern edge, assuming a discontinuous 
jump in the gas density, which produced a best-fit density jump of 
2.6 ±0.2. We also smoothed the best-fit model by convolving the 
density jump with a Gaussian function to determine if a non-zero 
width produced a better fit to the surface brightness edge. A model 
with f7 = provided the lowest and we calculated a 95% upper 
limit on the cold front edge width of (7 = 2 kpc. 

This upper limit on the width can be compared with the 
Coulomb mean free path of the electrons and protons on both sides 
of the cold front. As in section |3.1| we estimate the mean free 
path for particles crossing from inside the subcluster core to the 
outside, Ain->out = 0.5 — 1 kpc, which is the main source of dif- 
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Figure 17. Surface brightness profile for the southern edge of the subcluster 
core with the best-fit projected discontinuity model shown overlaid (solid 
line). The dashed line corresponds to the best-fit projected discontinuity 
model smeared with a Gaussian of width a = 5 kpc and the dot-dashed line 
shows the model for the outer cluster emission. 



fusion across the edge. If Coulomb diffusion is not suppressed 
by any mechanism the southern edge should have a width of at 
least several times Ajn^Quf Therefore it is hkely that transport pro- 
cesses are highly suppressed here. As discussed in the introduc- 
tion, the gas motion around the subcluster core is likely to produce 
a preferentially tangential magnetic field and strongly restrict heat 
flux and diffusion across the front (eg. |Vikhlinin, Markevitch & 



Murray|2001bl|Narayan & Medvedev 2001 



Asai, Fukuda & Mat- 



sumoto'2004' '2005 ' '2007 'Lyutikov 2006;"Xiang et al.'2007[ |Dursr 
& Pfrommer,2008, see also ^Churazov & Inogamov_2004| l. 



4.3 Structure in the subcluster tail 

The NE and SW sides of the subcluster tail appear strikingly dif- 
ferent in this deeper Chandra observation. Even in the raw counts 
image (Fig. [TJ, the SW side of the tail exhibits a sharp, coherent 
edge, that stretches --^150 kpc in length. Fig. [2] (lower left) indi- 
cates that there may be a second fainter edge, 100 kpc in length, 
outside the main edge. It's unclear whether these are different or 
related structures as the subcluster tail is clearly three-dimensional 
and structures are seen in projection. The narrow SW edge of the 
subcluster tail appears similar to the sharp, straight edges of the 
subcluster bullet in the Bullet cluster which appear unaffected by 
turbulence ( |Markevitch|2006||Markevitch & Vikhlinin|20(J7| l. Tur- 
bulence may be suppressed along the SW edge of the core by a 
continuation of the magnetic layer that strongly restricts conduction 
across the leading edge of the cool core (section | 4.2[ eg. [Lyutikov] 
[20061 IXiang et al.|2007[ [Pursi & PfrommerpOOsV 

In comparison, the NW edge is broad (~ few x 10 kpc across), 
and poorly defined as it is disrupted by the shearing flow of the am- 
bient ICM. Small-scale instabilities are expected to grow rapidly 
on timescales much shorter than the cluster passage time and widen 
the interface. If magnetic draping is responsible for the stability of 
the SW edge of the subcluster core and tail, then it is initially sur- 



prising that the eastern edge of the core is not similarly stabilised. 
However, as discussed in section [2.3. 2| the orbit of the subcluster 
core is likely to be curved to the South, which will increase the rela- 
tive velocity of the ambient ICM along the NE edge of the core and 
reduce it along the SW edge. The ambient gas motion is responsi- 
ble for the accumulation, stretching and ordering of the magnetic 
fields along the subcluster interface and an alteration in the subclus- 
ter trajectory could also affect the build up of this magnetic draping 
layer. There is no significant variation in the thermal pressure along 
either side of the subcluster core (Fig.[T4| left) but there could be a 
significant variation in the flow velocity of the ambient ICM caused 
by the subcluster trajectory, which is affecting the stability of the 
eastern edge. 

Fig. [141 shows the variation in the ICM pressure and temper- 
ature across the subcluster tail. The thermal pressure peaks in the 
dense, cool subcluster core but is also high behind the AGN where 
the flow of the ambient ICM converges behind the core. Although 
there is some fluctuation in pressure through the subcluster tail, and 
the eastern edge of the core appears to be breaking up, there are no 
significant structures consistent with large-scale hydrodynamic in- 
stabilities (eg. Abell 3667, |Mazzotta, Fusco-Femiano & Vikhlininl 
|2002| l. The cool blobs of gas breaking off from the eastern edge of 
the cool core (Fig. \T5\ appear to thermalize rapidly with the sur- 
rounding hot ICM and there is a steadily increasing temperature 
gradient through the subcluster tail. There is some significant vari- 
ation in the temperature gradient along the northern edge of the 
subcluster tail (Fig.|14^ where hotter 6 — 6.5 keV gas appears to be 
impacting a region of cooler 5 keV material. However, it's not clear 
if this relatively cool region has been stripped from the subcluster 
core and become isolated from the ~ 8 keV ambient ICM or if it 
has a similar origin to the SW plume and is a remnant of the pri- 
mary cluster core. Either way, this is likely to be a transient feature 
that will thermalize with the surrounding ICM. 



5 SUMMARY 

The deep Chandra observations of Abell 2146 have revealed a 
host of new and complex substructure. We confirm the detection 
of both a bow shock with Mach number M = 2.3 ± 0.2 and the first 
known example of an upstream shock, which has a Mach number 
M = 1.6 ±0.1. The bow shock, located ahead of the cool subclus- 
ter core, can now be traced to over ~ 500 kpc in length and ap- 
pears significantly broader than the upstream shock. We find that 
the observed postshock electron temperature profile behind the bow 
shock is lower than predicted for instant shock-heating of the elec- 
trons and appears consistent with the timescale for Coulomb coUi- 
sional equilibration with the ions. The equivalent measurement for 
the upstream shock appears to support the instant shock-heating 
model, similar to the result for the Bullet cluster. However, this 
measurement was more difficult for the upstream shock as the 
Mach number was lower, the uncertainty on the temperature val- 
ues was greater and there is significant substructure generated by 
the passage of the subcluster core. We therefore place less weight 
on the result from the upstream shock. 

In addition to the deviations from thermal equilibrium, we 
considered the possibility of non-equilibrium ionization in the re- 
gion immediately behind each shock front. However, the relatively 
low spectral resolution, low photon count rates behind each shock 
front and large radial bins compared to the ionization equilibration 
timescale make this a difficult measurement with the existing data. 

Even in the raw counts image, the upstream shock appears re- 
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markably narrow and well-defined over ~ 440 kpc in length. We 
calculate a best-fit width for the shock of only 6^3 kpc (95% er- 
rors), which is significantly narrower than the estimated mean free 
path and suggests that there is a significant suppression of Coulomb 
diffusion across the shock front. The bow shock appears broader 
with a best-fit width of 12^^ kpc (95% errors), however the mea- 
sured width of the shock fronts can be affected by deformations in 
the front shape which smear the edge in projection. Both the bow 
and upstream shock fronts appear to be consistently narrow across 
their measured lengths. 

The deep Chandra observation has also revealed a cool, dense 
plvune of material extending ~ 170 kpc in length in a direction per- 
pendicular to the merger axis. This plume is likely to be the remnant 
of the primary cluster core which has been pushed forward and lat- 
erally (SW) by the impact with the subcluster core. There does not 
appear to be a symmetric feature extending to the NE which sug- 
gests that the subcluster did not collide head on with the primary 
cluster. If the subcluster passed to the northern side of the primary 
cluster, simulations suggest that disrupted material from the pri- 
mary core would be mainly ejected in the direction of the observed 
plvune. However, a large impact parameter seems unlikely given the 
significant disruption of the primary core and the symmetry of the 
two shock fronts about the merger axis. 

The surface brightness images of the subcluster core show a 
sharp, leading edge corresponding to an increase in the gas den- 
sity by a factor of 2.6 ± 0.2 and temperature by a factor of at least 
3.4+Q g. The gas in the core has a metalUcity of 0.8 ± 01 Z© and 
is multiphase with significant temperature components at 2 keV, 
1 keV and 0.5 kcV. We find an upper limit of 50Mr;;, yr^' on 
the mass deposition rate in the subcluster core and determine that 
~ 80% of this cooling occurs in a bright filament. The leading edge 
of the subcluster core is particularly narrow with a 95% upper limit 
on the width of only 2 kpc, compared to the mean free path of 
0.5 — 1 kpc, so Coulomb diffusion is significantly suppressed here. 
We suggest that the motion of the subcluster core through the am- 
bient ICM has produced a magnetic draping layer which strongly 
restricts conduction and diffusion across this edge. 

The ram pressure stripped subcluster tail has interestingly dif- 
ferent structure along its two sides, which are both being sheared by 
the surrounding medium. Cool gas is being stripped off the eastern 
edge of the core into a 45 kpc filament which has a steadily increas- 
ing temperature along its length. This material then breaks off the 
filament into small blobs of gas which are directed into the sub- 
cluster's wake. The disrupted NE edge is broad at ~ few x 10 kpc 
across. In contrast, the SW edge is narrow and well-defined over 
~ 150 kpc in length. The magnetic draping layer around the lead- 
ing edge of the core could be stabilising the SW edge against tur- 
bulent instabilities but this appears to be less effective along the 
NE edge. For this non-zero impact parameter merger, the trajectory 
of the subcluster is likely to be curved to the southern, which will 
increase the velocity of the ambient flow aroimd the NE edge and 
could cause the observed gas stripping. 
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